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In this communication, we analyze the case of 3 + 1 dimensional scalar field cosmologies in the 
presence, as well as in the absence of spatial curvature, in isotropic, as well as in anisotropic settings. 
Our results extend those of Hawkins and Lidsey [Phys. Rev. D 66, 023523 (2002)], by including the 
non-flat case. The Ermakov-Pinney methodology is developed in a general form, allowing through 
the converse results presented herein to use it as a tool for constructing new solutions to the original 
equations. As an example of this type a special blowup solution recently obtained in Christodoulakis 
et al. [gr-qc/0302120] is retrieved. Additional solutions of the 3+1 dimensional gravity coupled with 
the scalar field are also obtained. To illustrate the generality of the approach, we extend it to the 
anisotropic case of Bianchi types I and V and present some related open problems. 



In the past few years, the ekpyrotic scenario has been proposed as an alternative to the standard inflationary 
cosmology [1]. In this interpretation, the big bang is viewed as the collision of two domain walls or branes described 
by Einstein's equations coupled to a scalar field. The scalar in this case parametrizes the separation between the 
branes. Hence, in this dynamical setup, it is of interest to understand the coupling of scalar fields to gravity. 

This role has been investigated rather extensively in 2+1 dimensional setups; see e.g., the earlier works of [2-4]. 
More recently the 3+1 dimensional case has become of interest; see e.g., [5-8]. 

These studies have in fact motivated the present work. In particular, in a significant recent paper [5], Hawkins and 
Lidsey have proposed a connection between the so-called Ermakov-Pinney (EP) equation and fiat, 3+1 dimensional, 
scalar field cosmologies. This is of particular interest since the EP equation is a very special, linearizable, nonlinear 
ordinary differential equation (ODE) that can be solved exactly if the underlying linear Schrodinger equation can be 
solved. We note in passing that one of the authors of [5] went on to use this approach to illustrate analogies between 
scalar field cosmological models and the dynamics of moments of the wavefunction of Bosc-Einstein condensates [9] . 
This point further highlights the importance of this ODE that has been recurring in a variety of different areas as 
diverse as nonlinear optics [10], elasticity [11], quantum field theory [12] or molecular physics [13]. For a recent review 
of the EP equation and its applications, see e.g., [14]. 

One of the purposes of the present short communication is to demonstrate how the methodology developed in [5] 
can be generalized in the case of non-zero curvature. We will also explicitly provide a converse result according to 
which, given a solution of the EP equation, the solution of a corresponding Einstein equation, coupled to a scalar 
field, can be derived. The use of the converse result will be explicitly demonstrated in a special case, that of a linear 
Klein-Gordon equation for the scalar field. The special solution for this case provided recently in [6] will be retrieved. 
More general solutions will also be constructed. 

As our gravitational model, we will use the 3-dimensional Friedmann-Robertson- Walker (FRW) metric. While other 
models such as the Brans-Dicke metric, see e.g., [15], are also popular, the FRW metric is widely accepted in the early 
universe scenario that is of primary interest herein (see the references mentioned above). The relevant line element 
(in co-moving coordinates and with a "cosmological" time choice) will thus read: 



In the metric of Eq. (1), a is the scale factor while c describes the curvature of the spatial slice and can be normalized 
to the values —1,0, 1 in the hyperbolic, fiat and elliptic case respectively. 

However, the above setting seems to restrict our considerations to an isotropic scenario. To illustrate the generality 
of the EP reduction and the usefulness of the corresponding technology, we also examine in the present work, an 
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anisotropic case, namely the one of Bianchi types I and V. In the latter, spatially homogeneous, yet anisotropic 
geometries, the line element is given by: 

ds^ = -N{tfde + 'ya0{t)a^{x)a^{x)dx'dx^ (2) 

in the time gauge N{t) = \/j(t)', o-f are the invariant basis one- forms of the homegeneous surfaces of simultaneity 
St, while are the scale factors that constitute the (in principle) dynamical variables of the metric. We restrict 
ourselves in the present study to these anisotropic models, as these have been argued to have the desirable property 
of isotropizing at arbitrarily long times [16]. Additionally, type V is the simplest Bianchi model that admits velocities 
or tilted sources, hence it is natural to consider it as the simplest extension that would allow the universe to choose 
a reference frame at the exit from inflation (given that the de Sitter metric does not have a preferred frame). 

Our presentation will proceed as follows: in section 2 we will provide the general EP methodology for the FRW 
system in the presence of a scalar. In section 3, we will apply these results in the special case of a massive scalar. In 
section 4 we will give additional special solutions to these equations, while in section 5, we will apply the method to 
the anisotropic case of Bianchi types I and V. Finally, in section 6, we will summarize our findings and present our 
conclusions. 



II. 3+1 ISOTROPIC SCALAR FIELD WITH CURVATURE: THE EP REDUCTION 



We will follow the notation of [5] in what follows. In particular, in this setting, the Einstein equations of gravity 
and the Klein-Gordon equation for the scalar field can be respectively written as 



dV_ 



0. 



(3) 

(4) 



H = a/a represents the Hubble parameter. The first two terms in the bracketed expression of Eq. (3) denote the 
energy density of the scalar field with potential F(0). The last term is the density of matter for the barotropic fiuid 
with equation of state Pmat = {n — 3)pmat/^] D > 0, < n < 6, = Sir/ml,, where mp is the Planck mass. Finally, 
the dots will be used for differentiation with respect to the cosmological or coordinate time t. Notice that with respect 
to the corresponding equations (2)- (3) of [5], the former is augmented by the spatial curvature term (recall that c can 
take the values —1,0 and 1 depending on the curvature of the hypersurface t =const.) in the left hand side. 

Now, by using a further difi^erentiation of Eq. (3) and the substitution b = a"/^ (cf. with Eq. (8) of [5]), we obtain: 




2 r 



nD 
362 



(5) 



If we now define a new comoving time r such that f = b, then Eq. (5) becomes 



fb 
dr"^ 



+ 



1263 



+ 



26^' 



(6) 



where the chain rule has been used and the functions in eq. (5) have been assumed as 6 = 6(t(r)) and = 0(i(r)). 
This equation can be compared with Eqs. (lO)-(ll) of [5], with the difference being evident in the inclusion of the 
curvature term (the last term on the right hand side of Eq. (6)). 

This prompts us to (briefly) discuss the EP equation which naturally arises in this context not only in the flat case 
of c = examined in [5], but also in the case of n = 2. The latter is a remarkable example of a nonlinear yet integrable 
ordinary differential equation (ODE) of the form: 



F" + Q(r)y = A 



(7) 



The particularly appealing feature of this nonlinear ODE is that its general solution can be obtained, provided that 
one is able to solve the linear Sclirodinger (LS) problem Y" + Q{t)Y = 0. For details on the properties of the EP 
equation, the interested reader is referred to [5,14] and references therein. Here we just mention its basic superposition 
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principle property. Namely, if the linearly independent solutions of the LS equation are Yi{t) and Y2{t), then the 
most general possible solution of the EP equation is given by 

r (r) = {AY^ + BYi + 2CYiY2)'^^ (8) 
where A, B and C are constants connected through 

and the Wronskian W = YiY^ - Y2Y(. 

It is of particular interest to note that the result of Eq. (6) can be used conversely for constructing solutions of 
scalar field cosmologies, using the EP equation structure and solutions. The converse result can be proved in the 
following form: Given Q and A = — < 0, let F > be a solution of 

d^Y A nc 

-— + QY=—^ + — — ^ (10 
Define a new time coordinate t such that f = Y{T{t)) and a = y^/n^ ^^yi as a new field (j) satisfying: 



Q, (11) 



with Q ^0. Finally, define a potential: 



12 fdYV Y^ fd(j)V D 3c 



\ dr J 2 \ rfr 

Then, if we consider the triplet (a(r(i)), ^(r(t)), F (0)), the latter satisfies the Eqs. (3)-(4). 

Let us now give a number of specific examples, where this construction scheme can be used to obtain solutions to 
the 3 + 1 dimensional scalar field cosmology equations. 

III. APPLICATIONS 

A. Flat FRW Metric and Massless Scalar 

In the absence of matter, we can set D = A = 0, and consider for convenience the case of n = 2. For A = 0, 
and for the flat FRW metric e.g., for c = 0, Eq. (10) becomes 

d^Y 



dr^ 

,2 J,2 , 



QY = (13) 



Assuming the general case V{(j)) = m (j) /2, the particular scenario of a massless scalar yields V{(j)) = 0, hence Eq. 
(12) becomes: 

'^^^^'^^-l^O (14) 



rfr2 Y \dT J 
whose solution yields: 

Y{t) = At'/^; (15) 

hence from Eq. (13): 

3 1 
16 72 

This, in turn, through the definition of f yields: 



= 7^72- (16) 



and, thus, finally from Eq. (11): 



r{t) = ( A'/hi (17) 



which is the same solution as obtained in Eq. (8) of [6]. 
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B. A Generalization: Non Flat FRW Metric Coupled with Scalar 



To demonstrate the generality of the technique also in non-flat cases with c ^ 0, we examine the example of n = 2, 

K = 1 and D > c/3 (which implies that A + c < 0). 
In this case, Eq. (10) becomes: 

Motivated by cases in which we are able to solve the EP (or equivalently the underlying linear Schrodinger) equation 
exactly, we choose Q{t) = 3/(16t^) (cf. Eq. (16)). In this case, using Eq. (8), we can find the general solution to Eq. 
(19) as: 

Y{t) = (At^ + Bt^ + 2Ct) ' , (20) 



where AB — C"^ = 4(A+c). For convenience, we use our freedom of coefficients to set A = B = 0, hence C = 2^/\X + c\. 
We thus obtain: 

Y{r) = {2Cr)i; (21) 

From Eq. (11), we derive: 

0(r) = ^log(r). (22) 

From the definition of f, it can be seen that: 

r = ^{t-tof. (23) 



Finally, from Eq. (12), it foUows that: 



!/((/.) = - =5e-^^ (24) 



T 



where S = 9c/8 - D/{2c) + 3/2. 

This result generalizes the corresponding result of Eq. (31) of [5] in the case in which curvature is present. 



C. An Example of Quadratic Scalar Fields 

As another example, we will use a case in which we use the EP approach in a reverse engineering way. In particular, 
in Eq. (10), we will postulate the solution and we will obtain the potential that is compatible with this solution. We 
assume that matter is absent, hence D = A = 0, for n = 4 and c > 0. We now demand that Y{t) = 2Bt. Then 

^(-) = (25) 

From the definition of f, we obtain that 

r{t) = ^e^^*, (26) 
where A is (without loss of generality) a positive constant. Then from Eq. (11), we obtain: 



where a is an arbitrary constant, while the potential 

V{4>)^^-^+B\4>-af (28) 

is quadratic in 0. Notice that in this example, the restriction of n = 4 can be lifted: in particular substituting 
A — > A^/^ and B Bn/A, we can perform the same calculation for any power n, but the expression for V{(1)) is 
functionally the same as the one of Eq. (28) . 

We remark here that the solution found in this subsection was previously identified by means of a different (than 
the EP) approach in [17] (cf. Eqs. (24)-(28) therein). 
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D. An Example of Constant Scalar Fields 



In the same inverse procedure spirit, another convenient choice of Q{t) is Q{t) = (for n = 2). Then the EP 
equation has the straightforward general solution: 

y (r) = (^T^ + B + 2Ct) , (29) 

with 

AB-C^ = \ + c=-'^ + c = \. (30) 

Consequently, 

r{t) = ^ exp {a^IH) + ^ exp (-A^^H) - ^. (31) 
It can then be immediately seen that the scalar (f) and the potential V are constants, while the scale factor 



a{t) = a(0) cosh (^A^^H^ + y a{Of - j sinh (^A^/H^ (32) 

with 

It should be noted that this is a generalization of the solution of Eq. (24) of [5] in the spatially non-flat case. Similar 
solutions have been obtained through direct calculations (i.e., instead of the EP methodology) in [3]. 



IV. ANISOTROPIC GENERALIZATIONS: EP REDUCTION FOR BIANCHI TYPES I AND V 
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A. Bianchi Type I 

Since the anisotropic case was not previously discussed in the earlier work of [3-5] , we discuss it here in more detail. 
In this case, the dynamical scale factors and the one-forms in the metric of Eq. (2) are given by: 

7a/3(i) = 

We can define the tensor: Fju^ = — SttT/^v where G^i, = R^v — \g^,i,R is the Einstein tensor and T^^ = 
4>;ti4';u — ^9tiv{4''°'<f';a + m'^(t>^) is the energy momentum tensor. Then, the quadratic constraint is the equation 
Fq = 0, the kinematic equation is given by = 0, while the Klein-Gordon equation for the field is given by 
•^llu ~ "^^^ = oc T.^^" = 0. Notice additionally, that the two integrals of the motion, namely 7i = — F| = and 
h = Fl - F| = yield B{t) = A{t)e''*/^ and r(t) = A(t)e^*/^. 

Solving the Klein-Gordon equation for (l)"{t) and substituting the result into dtFg =0 (as well as solving Fq = for 
(f) and using the resulting expression in dtFQ = 0), one is led to a dynamical equation for the remaining scale factor 
A{t) in the form: 



kX KA'jt) XA'jt) AA'jtf <i>'{tf A"{t) _ 

4 + A{t) + A{t) + A{tf 2 A{t) ^^^> 



Using now: A{t) = F(t)2/" and a change of variable r = / Q{t')dt', we obtain: 
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Hence, a choice of time reparametrization according to: 



^ = ^ + X+——^ (36) 
(which leads to fl{t) = 0e('*+'**^*F(t)^^+"^/"', where ^ > is a constant of integration), results in the form: 

nr).Qnr).^^. c = r.^-S^. or) 

Eq. (37) is of the form of Eq. (6) and becomes an EP equation for the choice of n = 6, for k = —A. 

In the more general case, the problem becomes extremely complex as t depends on r through dr/dt = fl{t) and fl 

is itself a function of the solution. Typically, this problem will not be analytically tractable (if combined, it leads to 
an integro-differential equation for Y(t)). However, this can be used as an inverse problem: we can postulate T{t), 
derive from it fl{t) and Y{t) and use them in Eq. (37) to derive the form of (pit). 

A simple example of the above methodology can be given as follows: let us assume that r = eJ'^^^''*' /{k + A), then 
n{t) = e('=+^)*. Hence, choosing 9 = 1, for n = 6 (the EP-like case for Eq. (37)) Y{t) = Y{t) = 1. Then, using the 
Eq. (37)), we find that 



and hence (e.g., choosing the integration constant C = log(K + ^)^2(1^^^)^ ^(^) 
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B. Bianchi Type V 

In this case, the scale factors and one-forms are, in turn, given by: 

( B{tf \ 
7a/3(^) = T{tf , at{x) = 
\ A{tf ) 

In type V, there exist three integrals of the motion, namely: Ii = f / — = 0, /2 = — = Q and = Fq = 0, 
which yield: B{t) = A{t)e''*/'^ as well as: r{t) = A{t)e-'''/^. 

Once again following the same path as before and substituting into dtF^ = 0, we obtain a single ODE for A{t) of 
the form: 




-6 A{t) A'{t) - 24 A{tf A'{t) + ^'"^ " ' 



,5 , 96A'{tr 



A{t) 

12 A{t) A'{t) ^'{tf - 24 A'{t) A"{t) = (39) 

Through a similar motivation as in the previous subsection, we use the change of variables: A{t) = F(t)^/", alongside 
the time reparametrization t = J* fl{t')dt', with fl{t) = ^y(^+")/". This results in the final form: 

•• ^ / \ '"■'^^ 1 n 1 

Y{t} + QY{t} = —^Y(Ty + -LVn ~ 2ff2 y(T-)l+4/n ' 

where Q = n , which is again a generalized form of an EP equation. This statement is made in the sense that 
while there is no choice of the exponent that yields a mere inverse cubic dependence on Y(t) in the right hand side of 
Eq. (40), the typical scenario involves inverse power dependence in a manner similar to the EP equation. 

While the resulting Eq. (40) is not directly of the form of the EP equation, its dynamics can be controUably 
tuned to be close to those of the EP solutions. In particular, it is clear that the asymptotic behavior of the equation 
will be similar to that of Y{t) + QY{t) = — fp- y^^^i"'"+i2/n for small initial data, while it will be close to that of 
Y{t) + QY{t) = ~ 2^ y(^)i+4/n for the case of large initial data. Hence, the choice of n = 6 for small initial data and 
the one of n = 2 for large initial data yields EP behavior (this point has also been verified in numerical investigations 
of the equation not shown here). In the intermediate regime where there is competition between the two terms the 
EP results will not be immediately applicable and one should resort to numerical simulations of Eq. (40). 
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V. CONCLUSIONS 



In this short communication, we have generahzed the earher work of [5] and of [4] in the direction of obtaining 
sohitions to 3 + 1 dimensional cosmological models for the Fricdmann- Robertson- Walker case in a systematic way. The 
method re-casts the relevant ordinary differential equations into one of the Ermakov-Pinney type which is explicitly 
solvable. From the solutions of the resulting EP equation one can re-construct the solutions to the original cosmological 
model in a step-by-step inverse process that has been detailed. The present study generalizes that of [5] in that cases 
of non-zero curvature of the FRW metric can be also considered. It also extends the results of [4] in the 3-1-1 
dimensional context. The method can be used to derive a variety of solutions in the latter context, including ones 
considered previously in [2,3,6]. 

On the other hand, we have also illustrated the generality of our method, by exploring the potential of such 
Ermakov-Pinney reductions to systems with anisotropy. In particular, we have demonstrated that the case examples 
of Bianchi types I and V can be reduced to Ermakov-Pinney-like equations. These Bianchi types were chosen as 
prototypical examples where a velocity/tilt can be included and as examples of Bianchi types that may eventually 
isotropize. The EP reduction of these cases highlights the generality and usefulness of the procedure. However, by the 
same token and since the EP is the only one of these inverse power, nonlinear ordinary differential equations for which 
an explicit solution exists, it underlines the importance of understanding the behavior of such classes of equations. 
In particular, examining numerically their temporal evolution for a number of physically motivated cases would be a 
natural next step. Such studies are currently in progress and will be reported in a future work. 
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